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We propose a non-perturbative method to determine the mixing coefficients of As = 2 four-quark operators for 
the Wilson quark action using chiral Ward identities. The method is applied to calculate Bk in quenched QCD. 



1. Introduction 

An essential step in the calculation of Bk with 
the Wilson quark action is the resolution of the 
mixing problem of the As = 2 four-quark op- 
erators, which is made difficult by the chiral 
symmetry breaking effects of the Wilson term. 
An apparent deficiency of perturbation theory 
for this problem has been well known [Q, and 
most calculations have tried to resolve the mix- 
ing non-perturbatively with the aid of chiral per- 
turbation theory pi. This method, however, has 
not been successful, since it contains large sys- 
tematic uncertainties from higher order effects 
which survive even in the continuum limit. Re- 
cently the method of non-perturbative renormal- 
ization|| has yielded a As = 2 operator with a 
good chiral behavior]^. However, the underly- 
ing mechanism of improvement in this approach 
is not quite apparent. 

Our aim is to calculate Bk with a method 
which explicitly incorporates the chiral proper- 
ties of the Wilson action, and to examine whether 
the result is consistent with that using the Kogut- 
Susskind action. In this report we propose a non- 
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perturbative method to resolve the operator mix- 
ing problem based on chiral Ward identities(WI), 
and report first results of a calculation of Bk car- 
ried out on VPP500/80 at KEK. 

2. Formulation of the method 

Let us consider a set of weak operators in the 
continuum {Oi} which closes under chiral rota- 
tion S a Oi = icfjOj. The continuum operators 
are given by a linear combination of a set of lat- 
tice operators {O a }, Oi — J2 a ^iaO a . We choose 
the mixing coefficients Zi a such that the Green's 
functions of {Oi} with quarks in the external 
states satisfy the relevant chiral Ward identities 
to O(a). The identities can be derived in a stan- 
dard manner |^| and take the form given by 

-2 P z A (E x P a Worn U k i>{pk)) 

+<3<<5;(0)n*^(P*)> 

-i £*(Oi(0) iAPkW^iPi)) + 0(a) = 
with pk the momentum of external quark. 

The four-quark operator relevant for Bk may 
be schematically written as Ovv+aa = VV + 
AA with V = sj^d and A = s^^ 5 d. To- 
gether with Ova — VA, it forms a minimal 
set which closes under A 3 chiral rotation. The 
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mixing pattern of these operators takes the form 

Ovv+aa/2=Z vv+aa (O + ziOi H h Z4O4) 

and Ova = Zvv+aa • Z5O5 where the lattice op- 
erators in the Fierz eigenbasis are given by Oq = 
(VV + AA) /2, O x = (SS + TT + PP) /2, 2 = 
(SS - TT/3 + PP) /2, 3 = (VV - AA) /2 + 
(SS - PP), 4 = (VV - AA) /2 - (SS - PP) 
and 5 = VA. 

Let us take four external quarks with an equal 
momentum p 2 — /i 2 . Let Tvv+aa and Tva be 
the sum of Green's function on the left hand side 
of (1) with external quark legs amputated. Us- 
ing the projection operator Pi corresponding to 
Oj, we can write Tvv+aa — T5P5 and Tva = 
TqPq + T\ Pi + ■ ■ ■ + F4P4 . We then have six equa- 
tions for the five coefficients Zi, • • • , z$, 

TJZw+aa = c + c\zi + ■■■+ c\zr = 0(a) (2) 

for i = 0, • • • , 5. We may choose five equations 
to exactly vanish on the right hand side. In the 
present analysis our choice is i = 1, • • • , 5. The 
remaining overall factor Zvv+aa is determined 
by the non-perturbative renormalization method 
of ref. S. We convert final results for matrix el- 
ements into those of the MS scheme with naive 
dimensional regularization (NDR) in the contin- 
uum at the renormalization scale fi = 2GeV. 



Table 1 

Run parameters. 
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L A x T 


24 3 x 64 


32 ;i x 64 
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300 
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a" 1 


1.95(5) 


2.65(11) 


3.41(20) 


m s a/2 
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0.0144(17) 
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Figure 1. Mixing coefficients Z\, ■ ■ • , Z4 for K = 
0.15034 at fi = 6.3. Vertical line corresponds to 
p w 2GeV. 

nation jackknifc procedure for all measured quan- 
tities. 



3. Parameters of numerical simulation 

In Table 1 we summarize parameters of our 
simulations. The lattice spacing is estimated from 
m p . At each (3 we employ four values of the hop- 
ping parameter such that the physical point for 
the K meson can be interpolated from data. We 
take degenerate s and d quark masses, and esti- 
mate m s a/2 from mx/nip — 0.648. 

For calculating Green's functions in Ward iden- 
tities quark propagators are solved in the Landau 
gauge for point source at the origin with the pe- 
riodic boundary condition. We extract Bk from 
a fit of plateau of the ratio of K°-K° Green's 
function of Ovv+aa divided by the vacuum sat- 
uration of the AA operator. For this calculation 
quark propagators are solved without gauge fix- 
ing for wall source at the edges of lattice using 
the Dirichlet boundary condition in the time di- 
rection. Errors are estimated by the single elimi- 



4. Results for Bk 

In Fig. 1 a representative result for the mix- 
ing coefficients is plotted as a function of exter- 
nal quark momenta. Data show only a weak scale 
dependence in the range 0.2<p 2 et 2 <1.0. We take 
values of coefficients at p « 2GeV, which falls 
within this range for our runs at j3 = 5.9 — 6.3, in 
the following analysis. 

We note that non-zero values for z 2 contrasts 
with the one-loop perturbative result z 2 = 0. 
Other coefficients agree in sign, albeit larger in 
magnitude, and approach perturbative values as 
/? increases. 

We show in Fig. 2 the a dependence of the ra- 
tio (K°\dvv+AA\K°)/\(0\P\K )\ 2 extrapolated 
to vriq = 0, which measures the contribution 
of chiral symmetry breaking terms. Results are 
plotted both for our method (WI) and with 
tadpole-improved one-loop perturbation theory 
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Figure 2. Effect of chiral symmetry breaking for 
Ward identity (WI) and perturbative (PT) meth- 
ods. 

(PT). The scalar density P in the denominator 
is perturbatively renormalized for both cases. A 
significant improvement achieved with the use of 
Ward identities is clearly seen, the ratio becom- 
ing consistent with zero even at lattice spacing as 
large as m p a ~ 0.2 — 0.3. 

Within the one-loop resolution of operator mix- 
ing chiral breaking effects are expected to appear 
as terms of 0(g A ) and 0(a). A roughly linear 
behavior of the PT values is consistent with the 
presence of the O(a) term. Also they linearly ex- 
trapolate to zero within errors at a=0. This sug- 
gests that 0(g 4 ) terms left out in the one- loop 
treatment are actually small. 

Our results for Bx(NDR, 2GeV) are summa- 
rized in Fig. 3. The WI method gives reasonable 
values even at a finite lattice spacing. Errors, 
however, are large, and a continuum extrapola- 
tion is difficult at this stage. 

In order to reduce statistical errors at each 
j3, we employ an alternative method (WI[VS]) in 
which the denominator of the ratio for extract- 
ing Bk is estimated from the vacuum satura- 
tion of Ovv+aa constructed by the WI method. 
While this method gives results different from 
those of WI at a^O, the discrepancy is ex- 
pected to vanish as a— >0. A linear extrapo- 
lation in a yields B K (NDR, 2GeV) = 0.59(8). 
This is consistent with a recent JLQCD result for 
the Kogut-Susskind action, ^(NDR, 2GeV) = 
0.587(7) (17) § . 
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Figure 3. Bk as a function of a. 

The tadpole- improved one-loop results (PT), 
if extrapolated linearly in a, give _Bff(NDR, 
2GeV)=0. 59(10), which agrees well those ob- 
tained with the WI or WI[VS] method. Uncer- 
tainties associated with a large extrapolation has 
to be resolved for assessing the reliability of this 
approach, however. 

In conclusion our results for Bk demonstrates 
the effectiveness of the method of chiral Ward 
identities for constructing the As = 2 operator 
with the correct chiral property. This makes us 
hopeful to achieve the goal of a precision determi- 
nation of Bk with the Wilson quark action with 
further improvement of our simulations. 

REFERENCES 

1. See, e.g., C. Bernard and A. Soni, Nucl. Phys. 
B (Proc. Suppl.) 9 (1989) 155. 

2. M. B. Gavela et al, Nucl. Phys. B306 (1988) 
677; C. Bernard and A. Soni, Nucl. Phys. 
B (Proc. Suppl.) 42 (1995) 391; R. Gupta 
and T. Bhattacharya, Nucl. Phys. B (Proc. 
Suppl.) 47 (1996) 473. 

3. G. Martinelli et al, Nucl. Phys. B445 (1995) 
81. 

4. A. Donini et al, Phys. Lett. B360 (1995) 83; 
M. Crisafulli et al, Phys. Lett. B369 (1995) 
325; M. Talevi, this volume. 

5. M. Bochicchio et al, Nucl. Phys. B262 (1985) 
331. 

6. JLQCD Collaboration (presented by S. Aoki), 
this volume. 



